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I. INTRODUCTION 





In the Aé€ronautical Engineering curriculum graduate 
jevel aerodynamics course, AE-4501, thea students are exposed 
+o two computer <crograms. One of thes2?, prepared by the 
Douglas Aircraft Company, analyzes the potential flow arcund 
three-dimensicnal wings but is limited to incompressible 
flow [Ref. 1}. The other progran, prepared by Cebeci, 
calculates the friction drag for two dimensional 
incompressible flow over airfoils [Ref. 2]. A serious 
defect of these pregrams is that they are not 
fat e-OL—-the-art computer frogra os. The Douglas program 
does not consider th2 effects of compressibility and the 
bourdary layer program, in additicn to being restricted to 
incompressible flow, does not predict the laminar to 


turbulent transition point. 


Av EACKGRCUND 

In 1980 the Department of Aeronautics at the Naval 
Postgraduate School acquired an intricatS computer program 
recently develoved by the Boeing Commercial Airplane 
Company. This state-of-the-art program calculates 


three-dimensional transonic flow over wings and bodies in 





both the cuter-inviscid flow regicn governed by <he 
transonic rotential egquaticn and the thin layer in which th 
meest order, ccmpressible boundary layér equations aré 


assumed tc ke valid. 


vp 


Th Boeing yrogram as received was designed +o be 
executed on a CDC 6600 or a CYBER 175 computer and was 
written using CDC FORTRAN IV extended language. This thesis 
therefore was primarily cencerned with th2 conversion of the 
orogram to FORTRAN IV extended ccempatible with the Naval 
Postgraduate School's (NPS) IBM 3033 systen. The large 
modular orcgram was divided so that th? potential flow 
analysis 2crtien could be run separately. Simplified 


instructions for use of the program were also prepared. 


ie VISCOUS/INVISCID SYSTEF OF PROGRAMS 

The Viscous/Inviscid Wing System (VIWS) of programs 
caiculates three-dimensional transonic flow over wings and 
meng body ccmbinations including details of the laminar or 
Meeoulcnt flow in ths three-dimensional viscous boundary 
layer. The flew field is calculated in ‘two overlapping 
regions: an cuter inviscid flow region governed by «the 
transonic rfotential equation, and a thin boundary layer in 


which the first order, three-dimensional, compressible 


10 





boundarv laver equaticns are assumed to hold and in which 
“he effects of surface teat and mass transfer can be 
computed. A list of the VIWS of programs is presented in 


Table I. 


TABLE I 


Viscous/Inviscid Wing Syst2m of Programs 









Program Name Description 


FIC27 Jameson-Caughey inviscid, 
+ransonic wing code 





AYQ11IN Reads gecmetry & velocity 


Bees GOnNStructs coordinate 
systen 


VWIN ctential flcw boundary 
ever interface 


AY11AC1 Three-diménsional boundary 
layer pregram 


INTERP Bcundar ee pOce mela lL 
fliow See erface 


Graphics display programs 


The basic sequence of calculations used by the VIWS to 













S@eaein matched viscous and inviscid solutions consists of an 
iterative loot in which tke inviscid outer flow analysis and 
the beundary layer analysis are performed sequentially. The 
iterative sequence is ccntinued until either convergence 

(satisfactcry matching) is achieved, or «he naximum number 

of iterations specified by the user has been performed. The 
VIWS programing sequence is shown schematically in Fig. 1.1. 
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Final boundary 
layer calc? 


Generate boundary layer grid 
Ist interface program 


Calculate final boundary layer flow 
1. Boundary layer program—upper surf. 
2. Boundary layer program—lower surf. 









Calculate transonic potential flow 
FLO 27 
Generate boundary layer grid 
Tst interface program 


Calculate 3-D boundary layer flow 
1. Boundary layer program—upper surf. 
2. Boundary layer program—lower surf. 






Plotting programs 


Set up 6* for next iteration 
2nd interface program 

1. Interpolation 

2. 5° under-relaxation 
3. Convergence test 










6° Yes 
converged 
fe 


Loop 
=ITMA 


Figure 1.1. Viscous/Inviscid Interaction Procedure 


memes rco-entiail flow is calculated for the bare wing during 
Ome LIrSt iteration. In subsequent iterations, the effect 


of the boundary layer flcw on the outer inviscid flow is 





felt as a modification tothe wing shape tkrough the 
addition of the boundary layer displacemant thickness. 
Convergence is recognized, and the iterations are stopped, 
when the maximum change between the new and old displacement 
thickness, expressed as a fracticn of the maximun 
displacement thickness, is less than the convergence 
*«olerance chcsen Ly the user. 

The VIWS utilizes the Jameson-Caughey transonic inviscid 
peea program FLO2Z/, tO carry cut the potential flow 
analysis. The boundary layer analysis is performed by a 
finite difference boundary layer prediction progran 
developed by the Being Cemmerciai Airplane Company. The 
basic theorv Eehind the bcundary layer program is contained 
in [Ref. 3]. A detailed descripticn of th2 VIWS of programs 
Oexcluding the potential flow program FLO27) is contained 
in (Ref. 41. A basic guide to the use of the VIWS of 


programs is contained in [Ref. 5]. 
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TI. POTENTIAL FLOW PROGRAM FLO27 


Eecause of the extensive length and number of program 
modules in the VIWS, the Fotential Flow Program, FLO27, was 
Singled out for cenversion. It was anticipated that FLO27 
would be run separately at first and recombined with the 
other crogqram modules at scme later date when czhese modules 
were themselves translated for execution on the IBM 3033 


computer. 


Ave RE-PROGRAMING 

The Potential Flow Pregram, hereafter called FLO27, was 
received cn magnetic tape and loaded into =he IBM 3033 mass 
stcrage system using the Jcb Control Languége (JCL) routines 
presented in Arpendix A. The magnétic tape contained twenty 
fe) total files in which «he formet was 9 track, 1600 CPI, 
unlakled. The card image format for the sixteen (16) 
program files is 80 characters per record and the four (4) 
Output files contain 150 characters per record. The program 
and outout files on the criginal CDC tape ars listed in 
feo le Il. 

The FLO2Z7 pregram was converted to FORTRAN IV extended 


suitable for execution on the IBM computer using the NPS CDC 


4 





TABLE 


CDC Magnetic Tape Filés 


Ul VAP aehots FLO27 Potential Flow + 


z /3194 AY11IN Reads gecmetry velocity data 
constructs See ec Ssysten 

Ba7 378 VWIN Potential flow boundary-layer 
inter face 

4 4/6611 A411A01 BBEC TS Tomer Chew poundary 
layer program 

a 7 vod INTERP Boundar y-layer potential flew 
abla terface 


T3770 | course 

| 

Procedures files 
Boeing McLean computer progren 


OUT4IIL | Output from boundary-layer, 
lower surface 

OUT4S11U } Output from poundary-layer, 
upper surface 


to ITEM conversicn guide [fef. 6]. The first step taken 





consisted cf crogram compilation using the WATFIV compiler 


with its extended error messages. The listing which was 


15 





produced flagged all areas of the program which required 
revision. Program changes were accomplished utilizing «his 
WATFIV listing. Some of the mcre general and repetitive 


changes are listed in Table III. 


AE Ee ieee 


FLO27 Re-Programing Changes 


CDC Code IBM Code Change 


ipcables: ier Roe 
IWRIT, EE 


WRITE ( IWRIT,600) WRITE (6 ,600) 
READ(IREAF,500) READ (5,500) 


READ te WRITE 7 or Changed to READ 14, 
REWIN “5 WRIT ieee) ie REWIND. 14 


All of this type eliminated 







Eliminated from program 





Hga< 

oO ao) 

ror 

> re 

ea | 

. 8 
rs 
= 
wu 
KH 
ae 


















The most difficult change +o make occurred with the CDC 
Buffer IN or Buffer OUT statements which were used in the 
program to transfer portions of a three-dimensional array 
into and out cf main memory. The Buffer routines reduce tk2 
memory size required to execute the program. This statement 
tv7<ze cccurred inthe main program and severai of the 
subrcutines. 
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The change réeguired to translate this statement is presented 
below with the CDC code preceeding the IBM FORTRAN. 

Purr eR OUT T(N3, 1) (G(1,1, 1) ,G(NX,MY,1)) changed to 

WRITE (N3) ((G(I,J,1) ,I=1,MX) ,J=1,MY) and 

EDFFER IN (N1,1) (G(1,1,M),G (MX,MY,M) ) Charged <o 

READ(N1,ERR= ) ((G(I,d,M) ,1=1,MX) ,J=1,MY) 

The variable ERR was assigned the GO TO statement number 
of the UNIT statemen* immediately following the BUFFER IN 
mere of code. AS an example, i1f the UNIT statement 
follcwine the EUFFER IN ccde was - If (UNIT (N1).GT.0.) GO TO 
151, “then the number 151 was assigned to variable ERR 
follcwing the €qual Sign. All CDC UNIT statements were 
eliminated from tha FLO27 source code per Table II. 

In addition to the pregram changes required to run FLO27 
on the IBM ccmtutzer, several lines of codes were added to 
modify the output format to a more usable form. A 
subrcutine, VERTEC, which calls the Versatec plotter was 
also Baasa to enhance prcecgram usefulness. Mies selottang 
routine is user controlled through an input variable and is 
ex lained in the next section. The modified FLO27 pzrcgran 


source code is presented in Appendix E. 
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To facilitate program data entry several input variables 
which had reccmmended values were initialized t2 these 
values within the Main pregram and the subroutine GEOM. The 
initialized input variables and their vaiues are presented 


in Table IV. 


TABLE IV 


Initialized Input Variables 


[xsca |S 
[seas 
[rcons [Sit 
a a 
rT 
rs 
[sree] 
Se FC 
a a 
xd 
[siz 
—astxe df 
[stuns [oo 
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A complete description of each input variable in Table 


IV can be found on pages 19 through 23 of [Ref. 5]. 


Bs PROGRAM DESCRIPTION 

The FLC27 program is a computer codé written to analyze 
+he transcnic flow over awing alcne or a wing ona 
cylindrical fuselage. It uses a finite-volune formulation 
so sclve the exact potential flow equation in conservative 
ZOE . In the development of the Equations, the basic 
assumptions are; St caae PlOw eno we neaceOG WOGK <Lan¢cfer, 


isentropic flow, irrotational flow, ano body forces and a 


cerfect gas. The velocity vector in cartesian coordinates 
is 
a ” “n 
V= ui + vj + wk (21) 
where u, v and w are tke velocity components. The 


continuitv equation, assuming steady flow, is 


o (Pu) pp) ~ 0 (22) 
Next a velocity potential 1s introduced such ‘that the 


velocity components are calculated as che gradient of this 


potential. 


u = ,v=$ , w = (Zea) 


uy 





ip 


Meee the  4ntrodiuction of the velocity potential, th 


continuity equation 2.2 becomes 
2h eH) + (0H) + S36 e4) =o (2.4) 


Assuming no heat or work transfer, the energy equation can 


be written as 


“6 = 4 
pee Goal) ee = T ag So Pe 


T 11 
p ~ 2 ~ 


(2.5) 


The flow is assumed to be uniform in the far field. On the 
surface of the bedy, the normal velocity component is zero. 
The velocities and densities of the near field are 
normalized using the free stream velocity and density, thus 
Yeo = 1 and boo = ee Using the assumptions that the flow is 
isentropic and a perfect gas, the energy equation 2.5 can be 


Shown to ke 
1 


= ene a 
0 = 11+ ee M (-¥) i (2.6) 


With equaticns 2.5 and 2.6 there ar2 only two unknowns, 
D and fe 2 They can be solved, subject tc «he boundary 
condition of flcw tangency, using a finite volume technique. 


The basic numerical scheme for «he solution is the 
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constructicn cf a mesh frcem small volume elements (cubes) 
Sf 

which are packed around tke wing or wing body configuration. 

The cubes in the computational domain are separately mapped 


-» distorted cukées in the physical domain by independent 


transformations from local coordinates X, Y and Z to 
Cartesian coordinates x, y and z. The mesh points are the 
verzices (ccrners) of the mapped cubes. The velocity 


potential and dersity are calculated at each verctex in the 


mesh. The pressure distribution can then be calculated fron 
on (2.7) 


In the event that tke local flow velocity becomes 
sucersonic and shocks occur, these are handled in the usual 
Manner bv insuring that: 

1) The tangential’ velocity components are equal on each 
side of the shock. 

2) Continuity is maintained by keeping the product of 
YU, constant across the shock (where J. is the normal 
velocitv component). 

3) Discontinuous expansions (corresponding to an 


NYexpansicn shock") are excluded from the flow field. 
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The assumoticn cf isentrofic flow along with the existence 
ef shecks rresents a contradicticn which can only be 
resolved by limiting the flow to very weak shocks for which 
sntrepy and verticity generation may be ignored. Thus; 
solutions will be valid only for Subsonic free stream 
velocities. 

The main three-dimensional array containing the 
Serentzal function data is stored on disk, and special 
unformated incutyoutput statements are used to bring planes 
of data into central computer memory and to store updated 
planes of data back Qn thesdssk. Unech> construction of the 
computaticnal coordinate system, a Joukowski transformation 
is used to transform the cylindrical fuselage to 2a vertical 
slit and then a sheared parabolic transformation is used in 
planes containing the airfoil sections. A detailed 
mathematical formulation of the potential flow analysis is 
contained in [Ref. 7]. 


Nie Peegram Input 





The inewe £0 FOZ) conSists of variables which are 


read with an &F10.6 FORMAT. Each input card has a titl 


iD 


Card which precedes it. Tito tesGala cottains the incu 


variakle name and effectively labels the input data for easy 
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reference. The title for each input data is placed in the 


same column as the input data it labels. The title cards 
are read with a 20A4 FORMAT. All numerical input valu2s are 
real numbers. The following data deck, listed card by card, 
is the minimum input data required for a simple wing 
analysis. Each "card" can be interpreted as one line of 
data on your terminal. A complete sample data set is 
presented in Acrendix C. 
CARD 1 The Run Title (64 characters maximum) 
CARD 2 Title card for the input variables 
ENG, ENZ, cobb cH end fPLrorT 
CARD 3 
Geils. 1-10 FNX ~- Number cf computational cells in 
the chordwise direction for th? initial mesh. 
MAX = 160/2**n, wkere n = FMESH -1. (See Cols. 31-40 
for FMESH) 
Cols. 11-20 FNY - Number of computational cells in 
the nermal eee en PeCcm tne alrcfoll suifiace for 
the initial mash. 
MAX = 16/2**n, where n = FMESH -1. 
SOlS ee iss0INZ -—- Number of computational cslls in 
the srcanwise direction for the initial mesh. 


MAX = 22/2**n, where n = FMESH -1. 
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Cols. 21-40 FMESH - Determines the number of times a 
program qenerated computational mesh is refined. 
EMte = Sen yd. 0, 2.0 OF eo) fOr. coarse, nedium or 
fine mesh. If 3.0 is selected the prcegram will 
calculate flow over the wing for the coarse mesh 
then half the mesh size (medium), recalculate, then 
half the mesh again (fine) anddoa final potential 
Liew Calculation. Output parameters are printed for 
each mesh size fer which calculations were 
performed. 


Cols. 41-50 PELOT - Output flag 


Or.0 NOm@taieGicpUt Wie NOUtT Prin te==plot.of Cp 
1.0 = Normal output with printer-plot of Cp 
at each computational mesh point for 
2ach wing section. 
2.0 = Normal cutput with Versatec plots of 
Cp versus X/C for each wing section of 
the final mesh. 
CARD 4 Title card for the input variables 
FIT, COVO and P10 


CARD S5-M Cne card for each computational mesh. Total 


numker of cards equal to 4 = FMESH. 
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CARD 6 


CanbD 7 


Cols. 1-10 FIT - A parameter which fixes the 
maximum number of iterations the program will us¢ to 
converge the velccity pctential to a specified 
tolerance (COVO). This parameter must be repeated 
for each mesh refinement. 

Cols. 11-20 COVO - Velocity potential convergence 
Grit sria. This input variable is also entered for 
each selected mesh. A value of 0.000001 is 
reccmmended. 

Cols. 21-30 P10 - This parameter determines the 
suksonic point relaxation factor for the specified 
mesh size. A value of less than 2.9 must be entered 
for each designated mesh. Pecommeandsd values are: 
WeGOVECE Coarse, 1.3 for medium and 1.2 for the fine 
mesh. 

fiewie Ca=d fOr “the input “variables 


FEAGH, YA, AL and CDO 


Cols. 1-10 FMACH - Free stream Mach number 


Cols. 11-20 YA - Yaw angle in deaqrees 


Cols. 21-30 AL ~- Angle of attack in degrees 


US 





Sand 8&8 


CARD 9 


Cols. 31-40 CDO -=- Drag ccefficient due to skin 


fr Geeen < Unless known, an estimated value of 0.01 
is reccmmended. 
Title card for the input variables 


ZOOM Noe oR EDP, DIRED and FUS 


Cols. 1-10 ZSYM - The wing planform symmetry 
eriqger. 

0.0 = Yawed wing, has no spanwise? symmetry 

1.0 = Swept wing, has Sfanwise symmetry 
Cols. 11-20 FNS -This input variable z¢lls the 
oregram the total number of wing sections you have 
selected to define the wing half span. The number 
must bre at least three (3) but not more than eleven 
(11) sections. 
Cols. 21-30 SWEEF - Leading edge sweep angle in 
deqrees. 
Cols. 31-40 DIHED - Dihedral angle in degrees. See 

ide 25.1. 

Cols. 41-50 FUS - Input the fuselage tadius. Enter 


0.90 for a wing-alcne case. 
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dihedral 


é uw oa 


/ eo, 


/ oS. soe 


anhedral (neg. angle} 


Figure 2.1. Dihedral Angle 


Data input cards from 10 through 15 are used for 

defining wirg rplanforms and section géometrics. For the 
east a secticn, all data cards from card 10 through card 
15 must be used. For the second and subsequent sections 
there is an cption for skipping the wing section defining 
data (cards 12 through 15) and copying the data from that of 
Smee ckrevicus section. This option is controlled by «he input 
Variable FSEC. If this cption is not used, data cards fron 
10 through 15 must be repeated for each wing defining 
Beet ion. The number of wing sections which are defined is 
insut with the variable FNS. Remember, up to 11 sections may 
be defined, anda minimum of 3 sections must be defined. 
All wing planform and section defining geometrics must be in 
consistent units. Wing planform and section defining 


,Uantities are tresented in Fig. 2.2. 
PA | 





XL 





YLbe-~chord 
Y Fi 
X 


Figure 2.2. Wing Defining Geometry 


CARD 10 itle card for the input variables 


Zo gee, VG, CHORD, TLHECK, AT end FS EC 
CARD 11 
Cols. 1-10 ZS - The section spanwise coordinates 
(Szar co atmerne centerline and werk 
out bo ard) 
Somes 1 i204 —. Secri0hn jiead:ng sedge xX coordinates 


Cols. 21-30 YL - Section leading edgs Y coordinat 


iD 


eoussso0-40 “CHORD = Section chord length 
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Seep 12 


Gelssmeti-5>0 EFHECK =" The thickness scaling factor 


can be used to scale all Y coordinates of the wing 
SecCcien. Thus percent thickness and camber are 
increased (or decreased) accordingly. Use 1.0 if 
no scaling is desired. 

Cols. 51-60 AT - Th2 twist angle of seach section 
(gecmetric twist) measured from the X axis to the 
Chord tone. A positive twist angle reduces the 
section angle of attack and gives "washout". Use 
0.0 f£c= no  twast . 

Cols. 61-70 PSEC - This is a flag which determines 
whether or not the pregram reads wing section 
defining geometry from a previous wing section or 
from new defining geometry. For the first section 
defined you must set FSEC tc 1.0. Following the 
first sectior, if you define new section geometry 
then use FSEC = 1.0. If you want the program to 
read the section geometry defined from the previous 
séecticn then set FSEC = 0.0. 

Title card for the input variable 


FN 


Ph) 





CARD 13 


CARD 14 


Cells. 1-10 FN - This variable contains the number 
ef pcints which define the upper and lower surface 
of the section. A maximum of 161 points may be 
used. 

Prelewcatd fOr tke 19put Variables 


XEF(I) and YP (Tt) 


CARDS 15-1 tec 15-N Total number of cards ¢quals N, 


where N = integer part of (FN+2)/3. 
The X and Y coordinates at each point are entered 


in pairs, three points to a card. (See Appendix C 


Celee l= lO exe (2) = X coordinet>s of the wing 
section point 
ii-2O wip (t)e= Y “coordinate of the wing 


secticn point 


21-30 defining X coordinate for next 
por ae 

371-40 demining Y coordinate for next 
point 

41-50 defining X coordinate for 


following point 
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5 1-6 0 d2opotiengeay COoEdanate for 
FOLMoOmw Ing POLDE 

The X and Y coordinates of the wing section defining 
coints must be entered starting with the upper surface 
railing edge peint and preceeding along the upper surface 
to the leading edge, and returning along the lower surface 
to the lower surface trailing edge point. Mey 2S very 
im.ortent to define the section leading edge with a larae¢ 
number of closely spaced points. Suggest at least 0.05 
snacing or less ketween X coordinate values from 0.1 X/C to 
the leading edge, X/C = 0.0. Each X and Y coordinat2 point 
is nermalized using the chord length for that section. 


Bereeicn defining geometrics are illustrated in Fig. 2.3. 


start 
Xp, Yp 


—_ 


— chord — 


Figure 2.3. Section Defining Geometry 
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@ynp 10 Title card centaining the words in Cols. 1-80 


END OF CALCULATICNS 
CARD 17 Title card for the input variable 

FNX 
CARD 18 

Cols. 1-10 FNX - This variable indicates the end 

of a set of calculations and must be set equal to 

0.0. Its puzpose is to indicat? that the progran 

has run td. completion. 

ee Prcgrapm Output 
Outout from the FLO27 program varies with the value 

of the inrut variable FPLCT. When FPLOT is set equal to 0.0 
a rormal output is produced. This normal output contains 
memecraer of occurance); refined input geometry data 
Sielud:ng trailing edge slope andangl2 calculations; 
-terative soluticn of the potential flow mesh; section 
characteristics and wing characteristics. The iterative 
ton, secticn and wing characteristic data are repeated 
for each mesh refinement requested. This, le shee on put 
Variakle FMESB is set equal to 3, these data are calculated 
and cutput three times. The last data in the normal output 


consists cf the non-dimensionalized chord (X/C) and vressure 


a2 





coefficient (Cp) data at each computational mesh point for 


each wing section calculated during the final mesh. A 
sample of the normal output data is presented in Appendix D 
and recresents the output data freom Appandix C input davté. 

Pemevataable FPLOT 1s set equal £9 1.0, the output 
data is increased considerably. This output contains the 
nozmal output pius a line printer-plot of theepressure 
coefficient at each computational mesh point for gach wing 
section. The line printer-plot is produced for each wing 
section of each mesh refinement. The length of the cutput 
data with FPLOT set equal to 1.0 can approach 6000 records 
d2cending on the number of mesh refinements requested. 
These vlots are of questionable value and, therefore, an 
alternate plotting program was developed. 

When the variable FPLOT is set equal to 2.0, the 
normal output data is preduced plus a Versatec plotting 
Subrcutine (VERTEC) is called. The subroutine outputs, vie 
eae Versatec slottsr, plots of Cp versus xX/C for each wind 
section of the final mesh calculaticns. This routine is 
Eee. y putting into plot form the Cp and X/C numerical data 
contained in the normal output. A sample of the Versatec 


Po. 1S cresented in Fig. 2.4. 
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Figure 2.8. Versatec Plot of Cp vs. X/C 
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iit TNRE ACTIVE TNPUT PROGRAM £rLOZ2/7IN 


The input data file required for ths FLO27 program is 
extensive. Frrors in input data FORMAT will cause progran 
errors at execution time. In crder to eliminats these 
errors and reduce the input data workload, @ computer 
+2rminal interactive program was written. This interactive 
program, called FLO27IN, is a user-frisndly way of creating 
an input data file fcr the potential flow wing analysis 
program FLO27. The FLO2Z7IN program source code is presented 
in Arprendix F. 

The interactive program, FLO27IN, when executed ask's 
qguesticns cf the user in crder tc censtruct and write to the 
user's "A" disk the required FLO27 input data file. The 
follcwing presents the step-by-step procedure for executing 
ame interactive program FLO27IN. 
aiae ei---Log on to any IBM 3033 interactive terminal 

with your user number and password. 
STEP #2---Once logged on and in the CMS operation mode 
type: 
CP LIK O247P 191 120 RR then hit ENTER. 
STEP #3---The werd PASSWORD will appear, Type and ENTER 


AERO 
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STEP ‘*#4---Type and ENTER 


ACC 120 D 
SreP #5---Type and ENTER 
LOAD FLOZ7IN (START 

The screen will display the header for the interactive 
orogran. Answer 2ach question presented. At the end of 
each guestion in parenthesis is the input data variable 
associated with that question and whether the input 
carameter is a real number (R) or an intager (I). Example: 
==> Fnter the free stream Mach number (FMACH): (R). FMACH 
is the invut data variable for the question. As you proceed 
meeough the FLC27IN progran, Opportunities to review and 
change input data will be presented. Should it become 
necessary to change your input data after completing the 
FLOZ7IN srogram, you can simply XEDIT the created data file. 

The FLO27IN program also inccrforat2s a library which 
Semcains the wing-section defining data fora number of 
Current wing shaces. A copy of this library is presented in 
Pe end:x B. This feature will be displayed during program 
execution by the use of a menu from which the user can 


select a pre-defined wing section or defin2 his own. 
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Ucon ceomele=ion of user inputs to the interactive 
vstoazram three additicnal data lines are automatically 
Boatten to the bcttom of the input file. They are: 

END OF CALCULATION 
FNX 
9.0 

mn additicn, Job Control Language (JCL) cards are 
mermeten +o the top and bottcm of the file. All JCL cards 
start with a // format. After FLO2Z7IN has run to completion 
ty-e and enter RELEASE 191 to release the aero disk which 
was linked while executing the PLO27IN program. The created 
data fil2 is written to the user's "A" disk with <filename> 
<file+vne> of FLO27 DATAIN. Additional changes can be made 


Simply by entering th XEDIT mode and editing the file. 
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Ven 7 se SATCH SY STEM EXECUTEON 


The votential flow prcegram FLO27 can be executed after 
the inout data file has been created. The batch processor 
is required fcr FLO27 execution because of the extensive CPU 
time needed to run the progran. While in the XEDIT meode, a 
standard JCB card must be added to the top of the FLO27 
DATAIN file orior to submission for job 2xecution. The JOB 
card has the form: 

//joktname JOB (nnnn,pppp) ,tident!,CLASS=J 
jobname = may contain up to 8 alphanumeric characters, 


the first of which must be alphabetic. 


rnnan your user number 


pppp = project number, assigned by professor 

megent* = contains the user's own identification 
information. A maximum of 20 characters may be 
contained within the single quotation marks. 

After adding the JOB card to your data file, you are 
rsady to execute the program. Type SUBMIT FLO27 DATAIN and 
Meess ENTER. Eatch runs are normally net worth waiting for. 
To inguire about the status of the job, enter INQ and the 
job name used on the JOB card or "logoff". If the system is 
busy and the maximum mesh size was selected, it may be 


s2veral hours before your job is run. 
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When the jok is run the output will be spooled to «the 
batch printer located next to the VM printar in the main 
Bomeuter building. The title at the top of the printout for 
batch jobs is the name entered on the JOB card. Lf We 
desired to have the program output data spooled directly to 
the terminal, it will be necessary to add one additional JCL 
card to *+he input data set. This card must be placed 
immediately following the JOB card and has tha forn: 

//*MAIN ORG=NPGVM1.nnnnP 
where nnnn = ycur user number 

Inserting this card in the input data will cause all 
program output to be spocled to the user's virtual reader 
mere it may be looked at, printed cr transferred to his "A" 
disk. To enquire as to whether information is in the reader 
Simply type RDR and hit enter, then follow the instructions 


on the screen. 
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Vic PROGEAM TEST RESULIS 


The FLC27 program was tested in three stages; (1) during 
the reprograming phase for conversion completeness, (2) 
after sucessful conversicn with suitable wing data for 


program accuracy and (3) during an AE-4501 class project. 


A. ACCEPTANCE TEST DATA 

To test and ensure that «he FLO27 program was converted 
to IEM compatikle Fortran without error, an acceptance test 
data set was used. The acceptance test input and output 
data was supplied with the criginal CDC program source code. 
After conversion of the FLO27 program *o Fortran suitable 
for the NFS IEM system, the acceptance test input data were 
run and «he output results compared to the output generated 
by tne CDC systen. 

Poth outout data sets were numerically exact when the 
FLO27 program was run in double precision on the IBM systen. 
If the oregram was run ir single precision, the numerical 
output values were exact tc the third decimal place. The 
difference in single precision accuracy occurs because the 
CDC system uses a 64 bit word length while the IBM system 


word length in single precision is only 32 bits. It was 
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decided that the IBM single precision accuracy was 


satisfactcry. 


B. CCMPARISCN WITH OTHER FROGRAMS 

The FLC27 program was also tested for accuracy by using 
the wing ovlanfcrm and section data from a2 NACA 572 wing. 
The data wer2 run on both the FLO27 program and the Douglas 
potential flow program [Ref. 1]. The data generated by both 
programs was ccmcvared to wind tunnel data for the NACA 572 
wing (Ref. 8}. The resuits are presented in Fig. 5.1 as 
plots cf lift ccefficient versus angl2-of-attack. The 
results shew that for the NACA 572 wing the FLO27 program 
more accurately predicts the wing lift coefficient than does 


the LCeuglas prcecgran. 


od AE-4501 CLASS PROJECT 

The final test phase was conducted by introducing the 
PLO27 program into the AE-4501 course as a class project. 
This was accomplished to determine student problems/comments 
concerning the data input program FLO27IN and to test an 
additional wing shape. The wing chosen for study was that 
of the A-7 airplane. The A-~-7 wing has a distinct leading 
edge metch at the approximate mid-Span. When the planfeornm 


geometry was run with the notch included the FLO27 program 
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Figure 5.171. Program Calculated and Wiad Tunnel Data 


Gee tO COmplstion but gave negative values for section and 


“Otal induced drag coefficient. The value for the lift 
co2fficient was low for the freestream Mach and 


angle-of-attack used. It was found that if the notch was 
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excluded from the wing gecmetry input data the program 


results were satisfactory both for induced drag and lift 
Seetficient. 

From the AE-4501 class experience it was determined that 
sharp wing clanform discontinuities cannot be handled by the 
program. If hcwever, the changes in shape are gradual, such 
as awing glovz, the program output appears to be 
Satisfactory. Such was the case with the acceptance test 
case data where the wing geometry was that of the F-8 


supercritical wing which incorporates 2 wing glove. 
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